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I. INTRODUCTION 

An uniform expansion of special functions is very useful representation of them which is used in many branches 
of science. It is well-known, for example, the Debay uniform expansion of the Bessel functions To obtain the 
uniform expansion one usually uses the complicate calculations which exploit an contour integral representation of 
function (see for example In this paper we use the WKB approach to obtain an uniform expansion for the 

Legendre functions. Previously, this question was analyzed by Thorne in Ref.Q by using different approach and in 
Ref. for particular case of the Legendre equation. We would like to note that this special case of calculations plays 
an important role in the so called functional methods which are at present the most powerful method (see pj). 

The organization of this article is as follows. First of all, in Sec|n]we reobtain the Debay formulas for the uniform 
expansion of Bessel function by using the WKB approach. In Sec lIIII we apply the same method to the Legendre 
functions and their derivative. The Appendix contains the list of the first four coefficients in manifest form. 

II. UNIFORM EXPANSION OF THE BESSEL FUNCTIONS 

In this section we reobtain the well-known 1] uniform asymptotic expansion for the Bessel functions of second kind 
I n (n\) and K n (n\) for large value of n. 

These functions obey to the following differential equation 

W" + \W = n 2 ( 1 + ^) W, (1) 



A V A 2 , 

where the prime is the derivative with respect A. 

Let us represent the solution of above equation as a series over small value of 1 /n: 



W = Ce nS - 1+So n ~ kuJ k (2) 



k=0 



with ljq = 1. Using this expression in Eq. Q we obtain the chain of equations 



s'-i = ^1 + ^, ( 3a ) 

S' = -^{S'^ + ^SU), (3b) 



2S'_ 1 { ~ x A 

Xu/„ V e A(A 2 -4) 



J k+i 



2 VVTTaV 8(1 + A 2 ) 5 /: 



w k , (3c) 



where s = ±1 and k = 0, 1, . . . . With new variable t = 1/ \/l + A 2 , the last equation may be rewritten in more simple 
form 

^fe+i = J (* 2 (1 - * 2 R)' + It 1 - 5i2 H> (3d) 
Z o 
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where the dot denotes the derivative with respect t. 

The first integral of the Eqs. © has the following form 

S-! = e(t?(A) + C_i), 



where 



So = -Jln(l + A 2 ) + C , (4) 

F E /•* 

"fc+i = -t 2 (l - i 2 )^ + - / (l-5t' 2 K(t')rft' + C fe+1 , 



A 



t?(A) = v^+A^ + ln (5) 
1 + vi + * 

To find the set of constants Cfe, fc = — 1,0, . . . we take the limit A — > oo in our expressions J2J, (0J and compare 
them with well-known asymptotic formulas |l( 



/n(nA )«^=L=e" A (l + 0(i)), K n (n\)*J^-e- n \l + 0(h). (6) 
V27rnA A V 2 ™A A 

Because of the next term of expansion is 0(1/ A) we have to set C k = for k > 1. Taking this into account we have 
the following expression for uniform expansion in the limit A — > oo: 

W w -L e e ™ A Ce enC - 1+Co . (7) 
V A 

Therefore, e = 1 corresponds to the uniform expansion of J„(nA) and e = — 1 to the K n {n\). For coincidence the 
expression Q with the asymptotic expansions 10 we have to set C_i = Co = and C = \j^phvn for e = 1, and 
C = ^JVphl for e = -1. 

Therefore, we arrive at the following well-known formulas for uniform expansion of the Besscl functions 



In{n\) = V^^^Tn-^), (8) 



Kn{nX) = V£ e_ "" (A) ^ (_nrfecjfe(<) ' 

fc=0 



fc=0 



where 



cu k+1 = \t 2 {\ - t 2 )uj k + I I (1 - 5t' 2 )u k (t')dt'. (9) 



2 

In order to find formulas for derivative of the Besscl functions we represent them in the form below 



-W' = C e «s-i+s J2 n- k Zj k . (10) 



n 

k=0 



Comparing the derivative of Eq. with respect A with above formula we obtain 

S—i = S—i, 

So = So + lnieS'^), 

C = eC, (11) 

ui k = uj k + ^t(t 2 - l)cj fe _i + et 2 (t 2 - l)d?fc_i- 

Therefore, with these expressions we arrive at the well-known formulas for uniform expansion of the derivative of 
the Bessel functions 

1 1 1 00 

-7>A) = -—-e^W^n-^W, (12) 
n \J2-nnt A 
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1 I 1 °° 



where 

1 



ZU fc = uj k + ^t(t 2 - l)w k -i + t 2 (t 2 - l)cj fc _i. (13) 



III. UNIFORM EXPANSION OF THE LEGENDRE FUNCTIONS 

In this section we employ the same approach for the Legendre functions. We consider the following equation 

(1 - x 2 )^" - 2.t*' - (n 2 ^ + y~~^ + 2^ * = 0, (14) 

which has appeared in context of quantum field theory in curved space-time 0.0 Here x G (—1, 1), ",7 and £ are 
real numbers, and the prime is the derivative with respect x. The particular case of this equation for £ = 1/8 has 
been considered in Ref. Q. 

The solutions of this equation are the Legendre functions first and second kind: 

m, q*[x\ (15) 

with index 

V = -\ + ^1-8^-4^2- (16) 

For £ = 1/8 these functions are called the cone functions Q . 

We assume n > and consider the following two independent solutions 



P n M = P» n [x], (17) 
(-1)* 
2 



i n M = ^(Q^N + QViN) 



1 -K\-A- 



2siii7r/x~ M 

They are real functions for arbitrary [i and obey the following Wronskian condition 

1 



To obtain the uniform expansion of functions (|17|l for large number n we represent the solution in the WKB form 
as below 

oc 

* = Ce nS - l{x)+So{x) n ~ k Mx) (18) 

k=0 

with ipo(%) = !■ We would like to note the difference of the uniform expansion in form 1)18(1 . which is over inverse 
degree of n, with that considered by Thorne in Ref.0. He obtained an expansion over inverse degree of (j, + 1/2 = 

y/l-8£ - 4n 2 7 2 /2. 

Substituting above expression in Eq. 114|) we obtain the chain of equations 



7 2 



S'_ j — ■ £ * 

So = -U^-Ar 2 \ 



(l-x 2 ) 2 l-rr 2 ' 
{ Six _ 2x 
2 ~ 1 -x 2 



^ = 2S 



4 



= -7T 



(1-X 2 )^ 

2 [VT+fF^) 



+ Vi^/c, fc> l, 



where e = ±1. 

The first integral of this chain has the following form 



S-i(x) 



7 arctan 



jx 



v/l + 7 2 (l-x 2 ) 



1, (l + a;)(l+7 2 (l-a:) + 
- In 



; (l-x 2 )) 



2 (1 - x){l + 7 2 (1 + a:) + y/l + 7 2 (1 - a; 2 )) 



S (x) = --\n(l +1 2 (l-x 2 )) 



ipk+i(x) = C k+1 (e) - - 



1 - x 2 



2 x /l + 7 2 (l-a; 2 ) 



+ e 





' 7 2 [ 




V 8 . 



2 - x 



,2 



5.7,' 



/2 



(l + 7 2 (l-a;' 2 )) 1 /2 



(1 + 7 2 (1 - x' 2 )f/ 2 (1 + 7 2 (1 - a;' 2 )) 5 / 2 
ipk(x')dx'. 



We have already set the constant Co = 0. This leads to redefinition the constant C, only. 
The formulas look simpler in terms of new variable 



v /l+ 7 2 (l-a; 2 ) 

instead of x. This quantity obeys to inequality: \v\ < |x| < 1. In terms of this variable we have 
„ . . f 1 , 1 — v 

o-i{v) = £<— — m- — ; h 7 arctan 7t) + G_i 



S (v) = 



2 1 + v 



1, 1+ 7 V 
— In — — . 

4 1 + 7 2 



/ t\ r £ (1-^ 2 )(1 + 7 V) ; f s 
Vk+i{v) = C k +i-- (i + 7 2) ^ v > 



8(1+ 7 2 ) Jo 



' / -^^ 2 + ^ 1 + ^- 1 ) 72( i + 7 v 2 ) 



(20) 



(21) 

(22a) 
(22b) 

(22c) 



In above formulas the dot denotes the derivative with respect new variable v. 

In order to find constants Ck we have to compare our formulas with exact expressions for the Legendre functions 
at a fixed point. For this reason we take the limit x — > 1 in our formulas 



« C 



1-x 



-en/2 



exp 



ne ( C_i — — ln(7 2 + 1) + 7arctan7 



and compare them with well-known expressions |3( for the Legendre functions at point x = 1: 



P n M=p- n [x] « - 



1 /1-x 



a/2 



(23) 



(24) 



9pN = ^(QpN + QViW) 



(re - 1)! (\- x 



-n/2 



Therefore, from Eqs. I|23(l . (|24|l we observe that e = — 1 corresponds top™[x] with C = 1/n!, and £ = +1 corresponds 
to q™[x] with C = (re - l)!/2, and 



C-x = — ln(l + 7 2 ) — 7 arctan 7 



(25) 
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for both signs of e. Furthermore, the coefficients ipk{v) must obey the following condition 

^fe(l) = 0. 



(26) 



Taking into account above formulas we arrive at the following expression for uniform expansion of the Legendre's 
functions 



f 



l + 7 Vi 1/4 



1 + 7 2 



> nS -^J2 n ~ k Mv) 



2, ,2 



l + Y v 



1 + 7 2 



1/4 



k=0 



-nS—i(v) 



J2(-n)~ k Mv), 



k=0 



where 



i>k+l(v) 



2"'(l+t;)(l + 7 2 ) 
(1 -w 2 )(l+7 2 i; 2 ) 



— 7 [arctan jv — arctan 7] 



7 



2(1 + 7 2 ) 
2 



(27a) 
(27b) 

(27c) 
(27d) 



8(1 + 7 2 ) 



Taking into account the same procedure as we used above for the derivative of the Bessel functions we obtain the 
following formulas for uniform expansion of the derivative of functions and 



1 d nr l 

--r-PM 

n ax p 



2, ,2 



1 +Yv 



I + 7 2 



3/4 



1 d nr ■ 



(re - 1) 



l + 7 Vl 3/4 l + 7 2 



1 + 7 2 



+ r 



k=0 

oc 

-ns^(v)j2(- n )-^ k ( v ) 



1-v 2 

k=0 

7 2 j)(l-D 2 ) (1 - V 2 )(l + ^ 2 V 2 ) 



2(1 + 7 2 ) 



I + 7 2 



(28a) 

(28b) 
(28c) 



The first four coefficients ipk and %p k are listed in Appendix. 

From the recurrent formula (|27d|) it is possible to find the value of the coefficients ip(v) for 7 — > 00. Indeed, 
comparing Eq. (|27d|) in the limit v — > and Eq. in the limit i — * 1 we obtain the following relation 



^ fc (0) = (-i) fc+ V(i). 



(29) 



Now wc represent formulas obtained in slightly different form which is close to expansion the Bessel functions. We 
set x — cose and 7 = A/ sine and use the asymptotic expansion for gamma function from Ref.[4j 



1 1 00 

Inn! = (re + — ) In re — re + — In 2ir + 



B2k 



1 



1 1 00 

= ("-2 )lnn - n+ 2 ln2 " + g2fc(2fc-l)re^ 



^ 2fc(2fc- 1) re 
^ 2Jfcf2J 



2fc-l 



where are the Bernoulli numbers. 
With these notations one has 



p^cosej 



27rre 



e nnJ2n- k i;+(v) 



k=0 



sin e 
Are 



<[cose] = J e -«^(-re)-V fc + M 



fc=0 



1 dpg[g] 

71 rfx |x=c 



2imt 



k=l 



sine 
Are 



An 



(30) 



• 2 ' 
sin e 
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FIG. 1: The plot of the relative errors Ap m = (pJJ — (Pm)™)/Pm versus of A for e = 0.1, £ = and n = 4. Here (pji)m is the 
uniform expansion of the Legendre function p™ up to degree m. 



n dx F= cose 



5^S^)¥ 5?? 

fc— 1 



where 



r/ = In 

f = 



A 



A 



\/l + A 2 + cos e sin e 
1 

, v — t cos e, 



sin e tan e 

arctan — arctan 



VTTa 2 " 
1 1 



l-8£ 



4n 2 A 2 



2 2V ' sin 2 e ' 



and the coefficients ip^ are found from relation 



oo 

E 

fe=0 



n = exp 



B 



21; 



2fc(2fc- l)n 



2k- 1 



Y^n k ip k (v) 



(31a) 
(31b) 

(31c) 
(31d) 



k=0 



by comparing the same degree of n in the left and right hand sides. 

The expressions (|3UH have the form similar to that for the Bessel functions expansion given by Eq. |JSJ. Furthermore, 
it is easy to see that in the limit e — > (argument of the Legendre functions tends to unit and lower index tends to 
infinity) the uniform expansion obtained is transformed to the uniform expansion of the Bessel functions below 



lim //VMcose] = i n I n (n\), 
lim ti- n q n a [cose] = C n K n (n\) 

as it should be according with well-known formulas Q 

T 1 

Jn{x), 

L ZJ 

lim z n Q 



(32) 



x 

cos — 

z 

X 

cos — 

z 



(33) 



\Y n (x), 



where x = inX and z — inX/e. In this limit the function 77 given by Eq. I|31a|) coincides with function ?/ J5| in the 
uniform expansion of Bessel functions: 



lim 77 = In — , 

Vl + A 2 + 1 



+ Vl + A 2 . 



(34) 



The numerical calculation of the relative errors Ap m = (p™ — (p^) m )/p^ and Aq m = (g™ — (g?J)m)/?u are plotted 
in Fig. 1 and Fig. 2 for different m — 0, 1, 2, 3 as function A, where (p!]) m and (g„)m are the uniform expansions of 
the Legendre functions p™ and g™ up to degree n _m . The difference is smaller the greater A. 
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In conclusion we would like to summarize the results. In this paper we obtain the uniform expansion for the 
Legendre functions [x] and g™ [x] given by Eq. (fTTj) for large indices n and ju = — \ + \ \J\ — 8£ — 4n 2 7 2 as a series 
over inverse degree on n. These expansions of the functions are given by Eq. (|27[1 and by Eq. I|28|) for their derivatives 
with respect of argument x. The coefficients of expansion may be found from recurrent chain of equations l|27djl and 
(|28c|) . The first four coefficients are listed in Appendix. 
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APPENDIX: MANIFEST FORM OF FIRST FOUR COEFFICIENTS. 

Below are the expressions for first four coefficients ipk and ipk m which we introduced for simplicity the following 
notations: 

1 x 

S = arctan[ 7 ] - arctan[ 7 w], £ = v= —======. (A.l) 

8 v /l + 7 2 (l-x 2 ) 



1p2 



V3 



7 7 2 + 1 
1 {6(\ 2 6£ 1 



2 7 2 + 3 w(7 2 - 1) 5w 3 7 2 



24 



24 



2 V7 



7 7 2 + 1 



27 2 + 3 v(j 2 -l) 5v 3 7 2 



24 



24 



C(-i + ^ 2 ) 

2( 7 2 + 1) 



(Y 



1 



4 7 4 + 84 7 2 - 63 w(7 2 - 1)(2 7 2 + 3) v 2 (9j 4 - 58 7 2 + 9) 



1152 



192 



128 



5vV(2 7 2 + 3) 77v 4 Y(Y - 1) 385vV 



192 



576 

6 I 7 J 2 U J 7 2 + 1 



1 



2 7 2 



1152 

V (Y 



3„,2 



24 



5u 7 
24 



8 A 

7 



( 7 2 + If 



47 4 + 84 7 2 -63 w(7 2 -l)(2 7 2 + 3) u 2 (9 7 4 - 58 7 2 + 9) 5w 3 7 2 (2 7 2 + 3) 
1152 + 192 + 128 576 

77w 4 7 2 (7 2 -l) 385w 6 7 4 l , C 



192 



1152 



r 



128 
27 2 + 1 v 2 
2 7 2 + T 



C 2 (i-^) 

2 7 2 ( 7 2 + 1) 



(A.2) 



8 



c 



+ 



( 7 2 + 1)2 
1 

(7 2 + I) 3 



2 7 2 + 7 v(3 7 2 - 11) u 2 (2 7 2 + 3) w 3 (44 7 2 - 29) 35w 5 7 2 



48 



16 



48 



48 



48 



1112 7 6 + 1116 7 4 - 918 7 2 + 5265 w(4 7 6 + 728 7 4 - 4323 7 2 + 711) 



414720 



9216 



w 2 (2 7 2 + 3)(9 7 4 - 58 7 2 + 9) w 3 (2005 7 6 - 37671 7 4 + 37566 7 2 - 2025) 



3072 



27648 



77w 4 ( 7 2 - 1)(2 7 2 + 3) 13w 5 7 2 (1053 7 4 - 3706 7 2 + 1053) 385w 6 7 4 (2 7 2 + 3) 



+ 



4608 

17017w 7 7 4 (7 2 - 1) 85085w 9 7 6 



15360 



27648 



9216 



82944 





= ipo 




= Ipl 






4>i 


= i>3 



1, (A.3) 
1 

12' 

1 , 1 

12 ^ 288' 

1 , 1 , 139 

— ^2 H w\ H • 

YT 288 51840 



i, 



+ 



7 7 2 + 1 
2 V7 



2 7 2 + 3 i;(3 7 2 + 1) 7w 3 7 2 



24 



24 



*C 1 

7 7 2 + 1 



2 7 2 + 3 w(3 7 2 + 1) 7v 3 7 2 



24 



24 



2( 7 2 + 1) 



1 



( 7 2 + l) 2 



6 \ 7 



+ 
+ 

^3 = 7 ( - 



+ 
+ 
+ 

+ 



(2 7 2 - 27)(2 7 2 - 3) i;(3 7 2 + l)(2 7 2 + 3) i; 2 (15 7 4 - 62 7 2 + 7) 



1152 



192 



128 



7w 3 7 2 (2 7 2 + 3) i; 4 7 2 (99 7 2 - 79) 455uV 



576 



i / 5C \ 3 i (K\ 2 1 



192 1152 

2 7 2 + 3 w(3 7 2 + l) 7w 3 7 2 



2 V 7 / 7 2 + 1 



24 



24 



7 



(2 7 2 -3)(2 7 2 -27) w(3 7 2 + l)(2 7 2 + 3) v 2 (15 7 4 - 62 7 2 + 7) 7w 3 7 2 (2 7 2 + 3) 



(7 2 + I) 2 



{ 1152 

w 4 7 2 (997 2 - 79) 455«V 



192 



1152 



192 

7 2 + 1 \ 2 7 2 2 



128 

, C 2 {l-v) 
^ 2 7 2 ( 7 2 + 1) 



+ 



576 



C 



( 7 2 + 1)2 
1 

(7 2 + I) 3 



2 7 2 - 1 w(3 7 2 - 7) u 2 (2 7 2 + 3) « 3 (447 2 - 25) 35w 5 7 2 
48 + 16 48 48 + 48 

III27 6 + 54367 4 + 1242 7 2 - 1215 v(12j 6 + 904 7 4 - 428I7 2 + 585) 



414720 



9216 



v 2 {2~f 2 + 3)(15 7 4 - 62 7 2 + 7) w 3 (28077 6 - 42897 7 4 + 37458 7 2 - 1863) 

3072 27648 
i; 4 (997 2 - 79)(2 7 2 + 3) llv 5 7 2 (152l7 4 - 4762 7 2 + 1241) 455«V(27 2 + 3) 



4608 



15360 



27648 



385w 7 7 4 (5l7 2 - 47) 95095v 9 y 



9216 



+ 



82944 



(A.4) 
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